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We report an ab initio study of the electric-field gradient tensor (EFG) at Cd impurities located at 
both inequivalent cationic sites in the semiconductor ln203. Calculations were performed with the 
FLAPW method, that allows us to treat the electronic structure of the doped system and the atomic 
relaxations introduced by the impurities in the host lattice in a fully self-consistent way. From our 
results for the EFG (in excellent agreement with the experiments), it is clear that the problem of the 
EFG at impurities in ln203 cannot be described by the point-charge model and antishielding factors. 

Introduction 

Perturbed-Angular Correlations (PAC) and other hyperfine interaction measurements are 
widely used experimental techniques that provide local information on the interaction of a 
probe-nucleus with the surrounding electronic charge distribution [1]. An interpretation of 
such measurements can lead to a detailed knowledge of structural, electronic and magnetic 
properties of a soHd [2]. One of the measured quantities, the quadrupole coupling constant Vg, 
is proportional to the mayor component of the electric-field gradient tensor (EFG) and 
the nuclear quadrupole moment Q. Due to the r ^dependence of the EFG from the charge 
sources, the EFG "felt" by a nucleus reflects sensitively the non-spherical electronic charge 
distribution around the nucleus. Therefore, the EFG is one of the most important clues for the 
study of the electronic structure in solids and the nature of the chemical bonding. 

Usually, the probe-atom is an impurity dopant, which is present in either trace or small 
amounts in the host crystal. For this reason, interpreting experimental results involves 
understanding chemical differences between the probe atom and the indigenous ion replaced 
by the impurity. For long time, approaches for a quantitative calculation of the EFG were 
often based on the point-charge model (PCM). Since such calculations do not account for any 
onsite polarization, Stemheimer antishielding factors were introduced, whereas chemical 
bonding was usually neglected. The value of such calculations was often quite limited. For an 
accurate calculation of the EFG, the electronic configuration of the host, perturbed by the 
presence of the impurity, has to be determined. This can be done in the frame of the density- 
functional theory. In this sense, in 1999 we began a systematic study using the Full-Potential 
Linearized- Augmented Plane Waves (FLAPW) method of the EFG in doped oxides [3-5]. 

Among the binary oxides, those that crystallize in the bixbyite structure were subject of 
several PAC investigations using ^^^Cd as probe [6]. Due to the high degree of ionicity of 
these oxides, it was assumed in the past that PCM was accurate enough to exactly predict the 
EFG in these compounds [7, 8]. We report in this work a FLAPW study of the bixbyite In203 
doped with Cd. Our calculations predict that Cd induce relaxations in the In203 host. The 
excellent agreement between our results for V33 and the asymmetry parameter rj and the 
experiments, demonstrate that this approach is suitable to calculate accurately structural and 
electronic properties of the system and to provide a deeper understanding of them. 

Method of calculation 

In203 crystallize in the cubic bixbyite structure (a=1.0117i A [9]). In this structure, the 
cations form a nearly cubic face-center lattice in which six out of the eight tetrahedral sites are 
occupied by oxygen atoms. The unit cell consists of eight such cubes, containing 32 In and 48 
O ions. Two nonequivalent cation sites, called C and D, both with 6-fold O coordination, 
characterize the structure. Site D is axially symmetric and can be described as an In 
surrounded by six O atoms at the comers of a distorted cube, leaving two comers of the same 



diagonal free. In site C the cube is more distorted 
and the O atoms leave free two comers on a face 
diagonal. The positions of all atoms are determined 
by four parameters: the coordinate u fixed the 
position of the C-type In (In(Q) ions, while the O 
positions are given by the parameters x, y and z [9]. 

In order to calculate from first-principles the 
EFG tensor at a Cd dilute impurity in In203 taking 
into account the structural and electronic effects 
introduced by the impurity in the host lattice, we 
considered the unit cell of In203 periodically 
repeated containing a single Cd-atom. Calculations 
were performed with the FLAPW [10] method as 
embodied in the WIEN97.10 code [11]. Exchange 
and correlation effects were treated using the generalized gradient (GGA) approximation [12]. 
In this method, the unit-cell is divided into non-overlapping spheres with radius i?, and an 
interstitial region. The atomic spheres radii used for Cd, In and O were 1.11, 1.05 and 0.85 A, 
respectively. We took for the parameter RKmax, which controls the size of the basis-set in 
these calculations, the value of 7. The correctness of the choice of these parameters was 
checked by performing calculations for RKmaj(=6 and 8. We also introduced local orbitals to 
include In-4p and 4d, 0-2s and Cd-4p and 4d orbitals. Integration in reciprocal space was 
performed using the tetrahedron method taking 26 k-points (Cd at site D) and 60 k-points (Cd 
at site Q in the first Brillouin zone. Once self-consistency of the potential was achieved, the 
forces on the atoms were obtained [13], the ions were displaced [14] and the new positions for 
Cd neighbors were obtained. The procedure was repeated 
until the forces on the ions were below a tolerance value. 
At each step, the Vu elements of the EFG tensor were 
obtained from the V2M components of the lattice harmonic- 
expansion of the self-consistent potential [15]. 

There is still an important point to be taken into account: 
the charge state of the impurity system, the single acceptor 
Cd^^ in In2^ Os^". In203 is a semiconductor with the oxygen 
/7-band filled. When a Cd-atom replaces a In-atom, the 
resulting system is metallic because of the lack of one 
electron necessary to fill the oxygen /7-band. Comparison 
of Fig. 2. a, 2.b. and 2.c. shows that the presence of Cd 
produces the appearance of Cd-d levels and impurity states 
at the Fermi level (£/). The question is if the real system we 
want to describe provides the lacking electron or not, and if 
this point is relevant for the calculation. Here we will 
present results assuming that the system provides the 
lacking electron (hence, the impurity level at Ef is in a 
charged state). To describe this situation we added one 
electron to the cell that we compensated with an 
homogeneous positive background to obtain a neutral cell. 

Results and discussion 

Initially, we performed calculations in undoped In203. Since the u-, x-, y-, and z-parameters 
are free internal parameters, which tumed out to strongly affect the EFG at both cationic sites, 
we refined the In(C) and O positions by force minimization. The results (m=-0.0336, x=0.3910. 
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Figure 1: Scheme of the oxygen distribution 
around hi atoms (black balls) in 11203. In the 
real structure the oxygen atoms are slightly 
displaced from the comers of the cube. The 
arrows denotes the V33 direction in each site. 
Results discussed in this work are referred to 
the indicated axes system (parallel to those 
defmed by the crystalline axes). 
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Figure 2: Density of states (DOS) for 
a) pure ln203; b) One Cd impurity 
replacing a D-type In atom in ln203. 
c) One Cd impurity replacing a C— 
type In atom in In203. In b) and c) 
the DOS correspond to the unrelaxed 
structures. Energies refer to the 
Fermi level {Ep). The arrows indicate 
impurity states in the valence band. 



3^0.1545, z=0.3813) are in very good agreement with both X-ray diffraction results reported 
in the literature [9]. We also calculated the EFG tensor at both cationic sites in pure In203. As 
can be seen in Table 1, the obtained result for the EFGs at sites C and D are in relatively good 
agreement with those experimentally determined by PAC at Cd impurities [16]. This similarity 
can be understand from the similar electronic structure of Cd and In atoms. 

After the study of undoped In203, we replace one D-type In atom (In(D)) by a Cd. As can 
be seen in Table 1, the substitution of a indigenous In(D) atom by a Cd impurity does not 
produce significant changes in the EFG. Again, this result can be understand from the fact that 
In and Cd present similar electronic structures. But, beside this similitude in the EFGs, the 
substitution of a In(D) atom by a Cd impurity induce not negligible forces on the Cd-ONN- In 
order to study the relaxation introduced by the impurity we considered the Onn displacements 
(assuming that structural relaxations preserve the point group symmetry of the cell in its initial 
configuration) until forces vanished. Displacement of atoms beyond the Onn are very small 
and none of the conclusions of this work are affected by it, thus we will not mention it in what 
follows. As can be seen in Table 1, the six Cd-ONN relax outwards along the Cd-0 directions. 
The displacement of each of these O atoms is of about 0.1 A (5% of the un-relaxed Cd-ONN 
distance). The magnitude of this relaxation was unexpected, and contradicts the assumption 
that Cd does not cause local lattice distortions in bixbyites [7, 8]. Concerning the EFG at Cd 
impurities located at cationic site D, the only effect of the relaxation is a small reduction of 
V33. After relaxation, the results for the EFG are in excellent agreement with the PAC results 
[16]. We want to mention that the EFG orientation in Table 1 correspond to the case of Cd in 
the bixbyites Er203 and H02O3 [8]. Due to the similarities between the structures, we assume 
that the EFG orientation in In203 should be nearby the same that in Er203 and H02O3. 

Table 1 : Mayor components, V33, of the EFG tensor at Cd site (in units of 10^^ V/m^), asymmetry parameter ?] and 
Vss orientation obtained in the FLAPW calculations (this work) compared with experimental PAC results at 300 
K. (i(Cd-ONN) are the distances in A from Cd to its Onn. hi the case of In203, (i(Cd-ONN) refers to the Iu-Onn 
distances. The experimental EFG tensor orientations correspond to Er203 and H02O3 (see text). 
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We performed a similar study for the case of a Cd impurity located at cationic site C of 
In203. Similar to the case of Cd at site D, the substitution of a In(Q by a Cd atom does not 
produce significant changes in the EFG tensor (see Table 1), but the presence of the Cd 
impurity induce a strong relaxation of its six Onn- The magnitude of this relaxation is different 
for each of the three pairs of distances Cd-ONN, but in average is similar to those produced at 
the site D (5%). As can be seen in Table 1, the structural relaxation changes the magnitude of 
V33, in the symmetry of the EFG tensor and a small reorientation of V33. The change in the 
symmetry of the EFG tensor reflects that the relaxation is non-isotropic. Finally, when the 
relaxed position of the Cd-ONN are taken into account, the agreement between our FLAPW 
results and the experimental ones [17] is excellent. 



It is important to compare our results with those predicted by PCM. The results of these 
calculations are presented in Table 2. As can be seen, the agreement between PCM and 
FLAPW predictions is very good for the case of site D. From these results, it could be stated 
that the ionic model can be used for the description of the EFG at cationic sites in In203. 
Nevertheless, the PCM predicts a negative sign of the EFG for site C, in contradiction with the 
FLAPW results. If the relaxed positions obtained by FLAPW are introduced in the PCM, the 
discrepancy in the sign of the EFG at site C is removed, but the agreement between the PCM, 
FLAPW predictions and experimental results for V33 at both sites and the symmetry at site C 
is not good. From this results, it is clear that the discrepancies between PCM and FLAPW (a 
method that take into account without the use of external parameters the electronic structure of 
the Cd impurity and the structural and electronic process that its presence induces in the In203 
host) can not be associated only to the structural relaxations. We have to conclude that the 
purely ionic model in combination with lattice summations and the Stemheimer factor fails in 
the description of the EFG tensor at cationic sites in In203. 



Table 2: V33, T], and V33 orientation obtained in our FLAPW calculations (relaxed structures), PCM, and PCM 
considering the structural relaxations predicted by FLAPW. 
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In conclusion, ab initio calculations using the FLAPW method successfully predict the 
experimental EFGs at Cd impurities in In203 and yielded quantitative information on the 
lattice relaxation around the impurity. From our results it is clear that the problem of the EFG 
at impurities in this semiconductor cannot be described by a point-charge model and 
antishielding factors. A proper theoretical description of electronic properties of Cd in 
bixbyites should consider self-consistently the electronic structure of Cd and the impurity- 
induced distortions in the hosts without external parameters or arbitrary suppositions. . 
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